How to measure diffusional decoherence in multimode Rubidum vapor memories'i 
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Diffusion is the main limitation of storage time in multispatial mode applications of warm atomic 
vapors. Surprisingly, there is scarce data for diffusion coefficients of Rubidium in buffer gases and 
available measurement methods are poorly suited to modern experiments. Here we present a sim- 
ple, efficient and direct method of measuring diffusion coefficients in gases by creating, storing and 
retrieving spatially-varying atomic coherence. Raman scattering provides a necessary interface to 
the atoms that allows for probing many spatial periodicities concurrently. We single out the diffu- 
sional component of depolarization and determine the diffusion coefficient. We found the normalized 
diffusion coefficients for Rubidium atoms in noble gases to be as follows: Neon 0.20 cm^/s. Krypton 
0.068 cm^/s. Xenon 0.057 cm^/s. We are the first to give experimental results for Xenon and we 
recommend this gas for multimode quantum storage applications. 

PACS numbers: 51.20.+d, 42.50.Hz, 42.50.Gy, 42.65.Dr 



I. INTRODUCTION 

In recent years warm atomic ensembles have been 
widely used in many applications in quantum optics and 
atomic physics. The most promising applications include 
quantum repeaters [1], quantum memories [2] and ultra- 
precise magnetometry [3]. They have also been shown 
to be an effective medium for four- wave mixing processes 
[4], Electromagnetically Induced Transparency (EIT) \u\ 
and slow light generation [6]. 

An advantage of using warm atomic gas is undoubtedly 
the simplicity of performing the experiment and obtain- 
ing large optical depths. However, it entails fast thermal 
motion of atoms. Atomic motion limits the efficiency of 
many modern systems that use multimode properties of 
atomic ensembles. At best this motion can be slowed 
down and made diffusive by addition of a suitable buffer 
gas. Thus, it is important to know the diffusion coeffi- 
cients for proper experiment design and interpretation. 

An important example of systems limited by diffusion 
is storing and retrieving transverse modes and images 
in Gradient Echo Memory (GEM) 0, B, in collective 
Raman scattering [9] or in Electromagnetically Induced 
Transparency (EIT) [10]. Typically diffusional motion 
of atoms in a buffer gas limits the storage time of im- 
ages down to tens of microseconds (ill, Jl2| or restricts 
the number of spatial modes retrieved J9|. In other ex- 
periments based on EIT diffusion broadens the EIT spec- 
trum [13]. Knowledge of the exact diffusion coefficient is 
required to design diffraction cancellation |l3 |. 

Based on the aforementioned examples we envisage 
that using multimode properties of atomic ensembles 
will become increasingly important in the nearest future. 
Given this fact, one of the most important defies will be 
to master diffusion processes so as to go beyond current 
limitations. Reading through the papers we found all of 
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them to have two features in common. In none of them 
the choice of the buffer gas type and pressure was well- 
founded, nor was the diffusion coefficient reliably mea- 
sured. There has been no attempt to optimize experi- 
mental conditions by controlling diffusion coefficients and 
knowledge of the value of the coefficient itself was rather 
doubtful. 

Most of the experiments with EIT and GEM listed 
above used Rubidium embedded in Neon buffer gas. This 
was probably due to the favorable ratio of diffusion co- 
efficient to spin depolarization cross section in this gas 
which made it standard in optical pumping and magne- 
tometry experiments. The values of diffusion coefficient 
of Rubidium atoms in Neon referenced up to date were 
measured 50 years ago with substantial effort and care 
UH, [iSI]. However, probably due to indirectness of meth- 
ods available those days, the diffusion values obtained 
in [15, 16] differ by more than 50%. Diffusion measure- 
ments for Neon had been conducted several times before 
the 1970s alongside measurements of spin relaxation cross 
sections [15-18]. The values given in these papers differ 
from each other up to 3 times. 

To our knowledge only one Australian group tried to 
use Krypton in the context of quantum memories in their 
GEM experiment [19]. This is a sensible choice when 
it comes to fighting diffusion, since it is almost 3 times 
slower in Krypton than in Neon at the same pressure. 
However, the Australian group does not quote any exper- 
imental measurements of diffusion. In their most recent 
paper [8] the authors measured the value of the diffu- 
sion coefficient using updated method of Franzen [20, zli|. 
They compared this measurement with the values they 
had inferred from the lifetime of the memory and with 
theoretical calculations, getting up to 500% differences. 
We take it as a motivation to perform measurements us- 
ing up to date method, all the more since we have been 
able to find only one paper from 1972 that gives the dif- 
fusion coefficient for Rubidium in Krypton [22] . 

We have been surprised by the scarceness and limited 
reliability of data on the diffusion of Rubidium in no- 



ble buffer gases. We have seen a considerable need for 
reliable results which, we believe, will help design novel 
experiments with the multimode use of atomic ensem- 
bles. Therefore we have decided to present a method 
for measuring diffusion coefficients that is easy to imple- 
ment, robust and up to date. We also believe that gases 
different from those that have been used so far can turn 
out to be even more effective. 

Given the latest applications of hyperpolarized Xenon, 
it would be useful to have reliable data on the diffusion of 
Rubidium in that gas as well [23]. Yet, it turns out that 
its diffusion coefficient has never been measured but it 
was only deduced from cross sections of velocity changing 
collisions. We have measured and list herein the values 
of the diffusion coefficients for Neon, which is the most 
popular, and Krypton we are the very first to give the 
value of this coefficient for Rubidium diffusion in Xenon. 

This paper is organized as follows. In Sec. II we sum- 
marize previous methods for measuring diffusion coeffi- 
cients based on examining spin relaxation and spot out 
their weaknesses. In Sec. Ill we introduce the principles 
of our method. Sec. IV describes in detail the experimen- 
tal implementation. Sec. V contains the experimental 
results. Finally, Sec. VI concludes the paper. 



II. METHODS OF MEASURING DIFFUSION 
COEFFICIENTS 

A. Previous methods 

We can affirm that problems with considerable uncer- 
tainty and scarceness of data for diffusion coefficients of 
alkali atoms in buffer gases result from complexity of the 
measurement methods themselves. There were designed 
primarily for the purpose of measuring the cross section 
of spin-exchange collisions. 

The pioneering method that involved measuring diffu- 
sion coefficients was proposed by Franzen in 1959 [15], 
and all of the subsequent contributions used the ini- 
tial conception unchanged or slightly improved. Franzen 
measured longitudinal relaxation time of pumped Ru- 
bidium atoms in Neon or Argon buffer gas. He applied 
a pump probe method through which he examined the 
change in absorption of the probe beam as a function 
of depolarization time and buffer gas pressure. In or- 
der to retrieve the diffusion coefficient from such mea- 
surements, Franzen made a few assumptions. He theo- 
rized that atoms depolarize owing to Rb-Rb collisions, 
Rb-buffer gas collisions, and fully depolarize in contact 
with walls. He also assumed the cell was a uniformly illu- 
minated cylinder. Consequently, he obtained the formula 
for average atom relaxation time as a function of buffer 
gas pressures. In the vacuum the relaxation was due to 
the atomic motion, while at high gas pressures it was 
mainly due to spin exchange collisions. Thus, by fitting 
dependence of the relaxation time on pressure with the- 
oretical formula, Franzen and his followers were able to 



find spin depolarization cross section and diffusion coeffi- 
cient. This procedure requires precise knowledge of tem- 
perature, Rb-Rb depolarization cross section, pressure 
and cell dimensions. The assumptions made at the be- 
ginning are virtually impossible to verify. Nevertheless, 
[15] had a profound significance for the future research 
dealing with diffusion of Rubidium atoms, since among 
many it presents both decent experimental results and 
sensible conjectures. 

Franzen's method [15], without further alterations, was 
used to measure diffusion coefficient in various buffer 
gases and conditions [l6|, |2^, |25| . The method was mod- 
ified later to measure not only longitudinal, but also 
transversal decay of polarization [13, [iM] • A certain in- 
novation introduced afterwards was the use of a spa- 
tially limited beam [20] and thereafter a laser beam [21] 
to probe polarization distribution in the cell volume by 
changing the position of the probe beam. 

By applying the Franzen's method one cannot mea- 
sure the diffusion coefficient for gases with a large spin 
depolarization cross section. This is why the diffusion 
coefficient for Krypton was measured only once [22] and 
there is no experimental data on the diffusion coefficient 
for Rubidium in Xenon. 

All of the papers mentioned above were based on the 
probing of spin relaxation evolution in the whole cell. 
They varied pulse duration, beam separation, pressure 
and the resulting diffusion coefficient was usually one of 
many fit parameters. Their results display large spread 
with discrepancies up to 300%. 

Recent attempts to calculate the diffusion coefficient 
from updated Franzen's method [2(/, zl] and GEM stor- 
age times also exhibit an appreciable spread [8]. Never- 
theless the Franzen and his followers are still cited and, as 
far as we know, there have been no measurements using 
significantly different techniques. 



B. Description of the new method 

Diffusion of alkali metals in buffer gases imposes a sig- 
nificant limitation on storage time in contemporary ex- 
periments with atomic ensembles. Attempts to repeat 
complex methods of Franzen and followers mentioned 
afore were not effective while no alternative methods have 
been presented. Thence, we propose a new, up-to-date 
and simple method for measuring diffusion coefficients. 

In our method we probe depolarization of fast spatially 
changing atomic coherence. Atomic coherence is blurred 
over time by diffusion as shown in Fig. [TJ Thus, we ex- 
pect high spatial periodicity components of the coherence 
to disappear faster than those varying slowly. By mea- 
suring decay rates for patterns of different periodicities 
we are able to retrieve the diffusion coefficient. Spatial 
periods of created atomic coherence patterns are much 
smaller than the cell size, therefore the influence of the 
walls can be neglected altogether. 

Patterns of atomic coherence are created through spon- 
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Figure 1. (Color online) (a) Planewave components of spa- 
tially dependent atomic coherence poi{r) act as diffraction 
gratings deflecting laser beam at different angles, (b) Com- 
ponents of high periodicity are blurred faster due to diffusion, 
therefore intensity of deflected light will decay faster for higher 
angle of deflection. During storage time r, pattern component 
of speciflc periodicity corresponding to wave vector K blurs 
with decay rate DK^, where D is the diffusion coefficient. 
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Figure 2. (Color online) (a) Atomic coherence poi(r) is cre- 
ated in spontaneous Stokes scattering, (b) Difference between 
wave vectors of laser fleld ki and scattering fleld component 
k2 is stored in the spatially dependent phase of the atomic co- 
herence, (c) Atomic coherence can be converted back to light 
in anti-Stokes scattering, (d) Atomic coherence phase is im- 
printed on anti-Stokes scattering as momentum conservation. 
Unlike the Stokes process, anti-Stokes scattering determinis- 
tically retrieves created coherence pattern on demand. 



taneous (Stokes) Raman scattering. Each pattern com- 
prises many plane-wave components with different peri- 
odicities. Those components decay at a different rates 
due to diffusional motion of the atoms. After a certain 
storage time a relative contribution of each plane-wave 
component to the pattern can be measured by driving 
anti-Stokes scattering. Then each plane-wave compo- 
nent acts as a diffraction grating deflecting driving laser 
beam as illustrated in Fig. [1] (a). By measuring the 
intensity of the Stokes scattering light as a function of 
deflection angle and time between pattern creation and 
readout, we can calculate the lifetime of different plane- 
wave components constituting atomic coherence. This 
represents a sizable amount of data containing informa- 
tion about homogeneous decay rate and the diffusion co- 
efficient. Values of those coefficients are easily retrieved 
through observing that only diffusion gives rise to an ex- 
cess in periodicity-dependent decay rate. Below we de- 
velop formulas describing the process. 



III. THEORY 



Diffusion of atomic coherence field 



and homogeneous depolarization with rate 70 0, [2^ : 



^Poi(r,t) 



-DV^Poi (r, t) - 70P01 (r, t) 



(1) 



Evolution of each plane wave component of initial am- 
plitude P{K) and wave vector K is described by a simple 
exponential decay at a rate 7(K) = 70 + DK'^. 

In general, the spatially dependent coherence poi (r) is 
a superposition of many plane waves of initial amplitudes 
/3(i^), each decaying at a specific rate: 



Poi(r,t) 



^'^P{K)e 



-DKH^iYi-v 



(2) 



K 



By measuring decay rates ^{K) as a function of wave 
vector K we can infer diffusional depolarization DK^ 
because it is the only contribution to the decay rate 7(i^) 
depending quadratically on K. 



B. Creation and probing of atomic coherence: 
Raman scattering 



Our method relies on creating, storing and probing 
spatially dependent coherence poi(r) between two long 
lived atomic levels |0) and |1). We observe the diffusion 
at the storage stage. As atoms move, coherence in a spe- 
cific point To will reshuffle its values with the neighboring 
points. Evolution of atomic coherence in the dark will be 
described by the equation of diffusion with coefficient D 



The spatially dependent coherence poi(r) between 
ground levels of Rubidium can be both created and 
probed by inducing Raman scattering. 

In Fig. [2](a) we illustrated the atomic levels involved 
in the Stokes scattering leading to creation of coherence 
poi(r) between levels |0) and |1). Upon spontaneous 
scattering of laser beam detuned from the |0) ^ |e) 



transition both the scattered hght and atomic coherence 
are created. We shah consider a simple case where the 
laser beam and the scattered light are planewaves with 
the wave vectors, ki and k2, respectively. The differ- 
ence between the laser field wave vector ki and the cre- 
ated photons wave vector k2 is accumulated in atoms as 
a spatial phase of atomic coherence. As illustrated in 
Fig. El^b) the atomic coherence created will be of a form 
poi(r) = /3e^^-^ where K = ki - ks. 

Such periodic atomic coherence poi(r) can work as a 
diffraction grating and defiect a laser beam. This is re- 
alized in the anti- Stokes scattering process presented in 
Fig. [2](c) in which laser beam detuned from the |1) ^ |e) 
transition is scattered at an angle. The spatial phase of 
atomic coherence poi(r) is imprinted back onto scattered 
photons as follows from Bragg condition. Provided driv- 
ing laser beam is the plane wave with a wave vector k'^^ the 
diffracted light wave vector will be k2 = k'^ — K as illus- 
trated in Fig. [2](d). Therefore by observing the intensity 
of the light lAsifi) scattered at an angle 0^ we register 
a signal which is proportional to the modulus square of 
the corresponding plane wave component of atomic co- 
herence |/3(i^)exp(-7(K)t)p, with e = K/k[ = K\/2iT. 



C. Averaging and retrieving diffusion coefficient 



The quality and reliability of the experimental data is di- 
rectly refiected in this last fit. 



D. Angular blurring at readout 

So far we have assumed the driving beam to be a 
plane wave. The finite size w of the driving laser beam 
in anti-Stokes scattering results in limited resolution in 
probing a wave vector space. The angular spread of the 
laser beam driving the readout will be transferred onto 
the angular distribution of the scattered light due to the 
momentum conservation even for a planewave atomic co- 
herence poi(r). Thus for any specific angle of observa- 
tion = K/k'i we detect scattered light originating from 
several distinct Fourier components of atomic coherence 
pattern. The contribution will come from the component 
of the wave vector K and its vicinity to the spread <j. We 
expect the spread a to be of the order of the inverse of 
the driving laser beam size 1/w. 

The result of the limited resolution is an overall in- 
crease of the observed decay rates 7^55. It can be esti- 
mated by convolving the storage time-dependent Fourier 
distribution of atomic coherence poi (r) with a Gaussian 
of a spread a, yielding: 



We create spatially varying atomic coherence poi (r) in 
the spontaneous Stokes scattering process, which popu- 
lates various planewave components randomly. Nonethe- 
less, the average modulus square of the excitation ampli- 
tude {\j3{K)\^) created right after the scattering is set by 
the driving pulse parameters and can be kept constant 
between measurement series. Therefore, we can calculate 
average intensity of the light scattered at a certain angle 
= K/k'i and for the given storage time r, incorporating 
Eq. 0: 

{Ias{0 = K/klr)) = r^{K){\/3{K)\')e-'^^^^^ (3) 

where r]{K) is an efficiency of readout, a priori depen- 
dent on K. The only factor that depends on the diffu- 
sion time r is the intensity decay factor e"^^^^^'^ , which 
provides direct information about the decay rate ^{K). 
Therefore for a given angle of observation = K/k[ we 
can infer 'y{K) from an exponential fit to a series of ex- 
perimental data taken for successive r. 

By repeating the decay fits for many K we can gather 
and then fit the expected functional dependence ^{K) = 
7o + DK^^ to obtain the diffusional coefficient D. In 
principle, we could use only two measurements of de- 
cay rates 7 corresponding to just two different direc- 
tions. Note that thanks to particular way of populating 
the spatially varying atomic coherence, we create and 
probe many wave vectors concurrently without altering 
the setup. Thence we obtain many points corresponding 
to a broad span of K vectors in a single measurement se- 
quence, which provides for a robust quadratic fit of 7(i^). 



-fohs{K) = {-fQ^2Dcj^)^DK^ 



(4) 



Note that the term quadratic in K in the above formula 
did not change, therefore the procedure to obtain the 
diffusion coefficient D remains unchanged. We only have 
to assure that the term 2Dcf^ is small compared to DK^ . 
This can be done by increasing the beam size w. 

In conclusion all measurements can be completed by 
varying only the one parameter - the diffusion time 
r equal to laser pulse separation while collecting scat- 
tered light on the camera. Data analysis require three 
straightforward steps: averaging, exponential and even- 
tually quadratic fit to obtain D. This makes the whole 
procedure relatively quick and simple to repeat. 



IV. EXPERIMENT 

In our case the levels |0) and |1) between which we 
create atomic coherence poi(^) are hyperfine split levels 
52S1/2, F=l and F=2 respectively. Spatially dependent 
coherence poi(^) is created and read with the use of Ra- 
man transitions between this levels. The scattered light 
is separated from much stronger stimulating lasers by a 
polarizer and additionally filtered out spectrally by Ru- 
bidium 85 filtering cell. The scattered light is registered 
by an Electron Multiplying CCD (Hamamatsu) camera 
sensor which is placed in the far field. The main part of 
the experimental setup is shown in Fig. [3l further details 
can be found in [9]. 
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Figure 3. (Color online) Central part of the experimental 
setup. Rubidium-87 atoms are mixed with the buffer gas in 
a glass cell inside a double magnetic shielding. Pump, write 
and read are the laser beams used to prepare atoms in the 
ground state |0), create and probe atomic coherence patterns 
respectively. Stokes and ant i- Stokes scattering are singled out 
on polarizing beam splitter PBS and observed in the far field 
on the electron multiplying EM CCD camera. 
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Figure 4. (Color online) Single operational sequence. First 
Rubidium atoms are pumped on D2 line to F=l ground level 
|0). Then spontaneous Stokes scattering is driven to create 
spatially varying atomic coherence. Next we store the coher- 
ence in the dark for time r and we let it fade due to diffusion. 
Finally we probe the coherence driving ant i- Stokes scattering 



The operational sequence is shown in Fig. |4j Laser 
pulses were formed with the use of acousto-optical mod- 
ulators. We began with optical pumping of Rubidium 
atoms to the ground state F=l. The pump laser operates 
in resonance on D2 line. Then we created random pat- 
terns of atomic coherence in the Stokes scattering process 
driven by 2-5 fis pulses of the write laser detuned from 
the F = l^F' = 1 Dl transition line by 1 GHz to the 
red. The Stokes scattering was recorded with a camera 
for diagnostic purposes. Next the atoms were left to dif- 
fuse for time r. Finally we used 5 /is long pulses of read 
laser detuned 1 GHz from the F = 2^F' = 2 Dl tran- 
sition line to the blue to probe blurred atomic coherence 
/)oi(r) and record anti-Stokes scattering intensity in the 
far field Ias{0). 

Beam diameters 1/e^ and powers for the Stokes and 
anti-Stokes drive laser were 5 mm, 4 mm and 16 mW, 7 
mW respectively. They were chosen in order to achieve 
both good resolution as discussed above Eq. (H]) and a 
sufficient scattering rate (9|. 

We repeated the create-store-read sequence multiple 
times and recorded random patterns of anti-Stokes scat- 
tering, changing the storage time r. 

We used four Rubidium 87 cells with different buffer 
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Figure 5. (Color online) Average number of photons per pixel 
{lAs(0,r)) in anti-Stokes scattering process as a function of 
the storage time r and angle of observation 6. Data taken for 
a cell with 1 torr Xenon. The decay is conspicuously faster 
for higher angles of scattering. The contour lines show 1200, 
800, 400, 50, 10 photons per pixel. 



gases. Those were: Neon - the most popular and the best 
explored, at the pressure of 2 torr. Krypton - which had 
been used in similar cases and we had some data con- 
cerning its diffusion coefficient, at 0.5 torr and 1 torr as 
in: [19] and [8] respectively. Xenon - though its diffusion 
coefficient had not been measured before, we presumed 
it to be the smallest (due to the specific weight of the 
gas). All cells were 10 cm long, 2.5 cm diameter cylin- 
ders made by Precision glassblowing. The longest times 
for anti-Stokes scattering observations were ca. r =50 
/is, which corresponded to an RMS atomic displacement 
of about 1 mm, far less than the cell size. 

The cells temperature was stabilized at ca. 70° C, 
which corresponded to an optical depth of 40 and a con- 



centration of rubidium atoms 



10 cm '^. Cells were 



heated with bifilar windings, but the heating current was 
interrupted for the time of impulse sequence inducing 
Raman scattering. The cells were placed inside a double 
magnetic shield. 

For each diffusion time r we averaged 500 images of 
anti-Stokes scattering each time obtaining smooth sym- 
metric profiles. We subtracted the averaged background. 
Given that the most important thing for us is the inten- 
sity as a function of azimuth angle 6>, we also carried out a 
polar averaging around the laser beam direction increas- 
ing signal to noise ratio. The result of the measurements 
was average scattering intensity {lAsifi^r)) as a function 
of angle 9 and storage time r. 



V. RESULTS 

In Fig. Owe present a typical map of the averaged an- 
gular profiles of anti-Stokes scattering in cell with Xenon 
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Figure 6. (Color online) Average number of photons per pixel 
{Ias{0^t)) as a function of storage time r with exponential 
decay fits observed for four different angles 0. Decays are 
visibly faster for higher scattering angles. Data corresponds 
to horizontal dashed lines in Fig. [5l 



vs free depolarization time (/^5'(6>,r)) . The results are 
given as a function of the scattering angle 6>, which is pro- 
portional to the wave vector of the corresponding Fourier 
components of atomic coherence patterns = K\/2t:. 

In Fig. [6] we plot an average number of photons as a 
function of storage time r at a few angles marked with 
dashed lines in Fig. [5l The decay rate is faster for higher 
angles of scattering. The error bars in Fig l6] correspond to 
Icr uncertainty and were calculated from the full statistics 
of camera counts. We fit exponential decays to the data 
taking the error bars into account. It is worth underlining 
that the data fits well to the curve at each scattering 
angle. As shown in Fig. El for high scattering angles the 
average signal {Ias{0, t)) is at the level of one photon per 
shot which increases the uncertainty of 7 considerably. 



and repeating measurements at different temperatures. 
The spread of diffusion coefficients obtained was about 
5-12%, depending on the cell. We attribute this spread 
to beam wander, laser power and frequency instability 
and the drift of the temperature of the cell during each 
measurement sequence. 



Buffer gas Pressure 


torr 


D cmVs 


70, 


obs kHz 


Ne 2 




91 ± 11 




38 


Kr 0.5 




136 ± 9 




71 


Kr 1 




57.5 ± 3 




28 


Xe 1 




52 ± 3 




33 



Table I. Fit parameters for measured cells. Temperature: 70° 
C . 



In Tab. |T] we summarize fitted values from the mea- 
sured data from charts in Fig. [3 Note, that for mea- 
surements in Krypton the ratio of obtained diffusion co- 
efficients is close to the nominal pressure values. 

To further verify the accuracy of our results, we car- 
ried out reference measurements for Krypton at 1 torr 
using write and read beams reduced ca. 3 times, so that 
their diameters 1/e^ were 1.6 mm i 1.4 mm respectively. 
This time the measurement was definitely less accurate 
due to the spread of wave vectors of the read beam and 
due to aberrations in the imaging system. The diffusion 
coefficient measured lay within the range of 40 cm^/s to 
65 cm^/s, which is consistent with other results . 

Finally let us note, that the observed decay rate at 
K = 0, 70, 06s summarized in Tab. |T] is dominated by 
excess contribution due to finite read beam size, 2Da'^ 
in Eq. (Hj). From respective collisional cross sections [16] 
we estimate ca. 1.5 kHz of decay rate due to Rb-Rb 
collisions and contributions less than 150 Hz from Rb - 
buffer gas collisions. 



A. Retrieving diffusion coefficients 

Having obtained an exponential fits for decays in all 
directions, we can analyze decays rates jobs{^) as a func- 
tion of the deflection angle 0. In Fig. [TJa) we present 
fitted decay rates for measurements in cells filled with 
1 torr of Xenon as well as with 2 torr of Neon . In Fig. 
[71(b) we give measured data for cells with Krypton at two 
pressure values: 0.5 torr and 1 torr. 

As expected, the data fits to the quadratic dependence 
7o6s(^) = 7o,o6s + D0'^{27t/X)'^. Now our la error bars 
correspond to respective confidence bound from the ex- 
ponential fit. Note that for high decay rate values the 
uncertainties are large because they correspond to small 
and noisy signals. 

The measurements described above were repeated sev- 
eral times in order to make sure the results were reliable 
and repeat able. We changed the amount of the light 
generated by altering pulse duration of the writing laser 



B. Normalized diffusion coefficients 

To conclude, we shall present diffusion coefficients Dq 
normalized to the standard conditions, i.e. at atmo- 
spheric pressure and at 0° C. They are related to the 
measured diffusion coefficients as follows: 



Do = D{ 



760 Torr^ 



(5) 



Where P stands for gas pressure at Tq = 0°C, T - tem- 
perature upon measurement. Diffusion coefficient will 
scale as T^/^ according to [27] 

The temperature is known with good accuracy and the 
main error in determining Dq results from inaccuracy of 
gas pressure in cells, specified by the producer not to 
be worse than 10%, and from the spread of measured 
D values. Normalized diffusion coefficients are listed in 
Tab. [n] together with data published previously. 



1.5 



0.5- 











/ 


/ 




— Ne2torr, D=91(4)cm^/s 
— -- Xe 1 torr, D=53(4) cm^/s 




[ 










. I -1-^ 

















1.5 



0.5- 











v 




— Kr 0.5 torr, D=136(6) cm% 
— -- Kr 1 .0 torr, D= 56(3) cm^/s 


/ 














/ 






^ 










y _^ 


.-^ 




f-— ■» — ^'T' 






■ 





0.2 



0.4 o.e 

e [mrad] 

(a) 



0.8 



0.2 



0.4 0.6 

e [mrad] 

(b) 



0.8 



Figure 7. (Color online) Decoherence rates ^{0) as a function of scattering angle corresponding to the atomic coherence wave 
vectors K — 27r^/A. We fit the data with quadratic dependence jobs{0) = 70, ofes + DO'^ (2ti / X^ . Quadratic term of fit gives 
diffusion coefficient D. (a) Data taken for 1 torr of Xenon and 2 torr of Neon (b) Data taken for Krypton of different pressures 
0.5 torr and 1 torr. 



Buffer gas i^o[cm^/s] - this paper 



Do [cm /s] - previous results 



Ne 


0.20=b0.02 


Kr 


0.068±0.006 


Xe 


0.057=b0.007 



0.11 [nj, 0.18 [18], 0.31 [15, 17], 0.48 [16| 

0.1 [22] 0.007 [8] 

No experimental data 



Table II. Measured diffusion coefficients of Rubidium atoms in noble buffer gases: Ne, Kr, Xe. Results are normalized to 0°C 
and 760 torr pressure. 



The results we obtained do not differ from those ob- 
tained beforehand, which - as we can see from Tab. HH - 
were characterized by a noticeable spread. 



VI. CONCLUSIONS 

We have demonstrated a novel method for measuring 
diffusion coefficients of atoms in buffer gasses tailored to 
atomic memory applications. The method is based on 
creation of spatially varying atomic coherence fields, let- 
ting them diffuse in the dark and probing them. The 
fields comprise various spatial periodicities evolving con- 
currently, created and probed with Raman scattering. 

A diffusion is observed in the dark and it is the only 
mechanism to produce periodicity dependency in the de- 
cay rate. With our approach, the evolution of the compo- 
nents of the atomic coherence with various periodicities 
can be observed in parallel and the periodicity depen- 
dency in the decay rate can be singled out easily. This 
leads to quite direct determination of the diffusion coef- 
ficient founded on a basic time and angle calibrations of 
the experimental setup. However, our method also re- 
lies on using lasers of a few MHz frequency stability and 
magnetic shielding. 

We have verified that the method is accurate and re- 
peat able. Results are based on multiple independent 
measurements for a number of various K vectors which 
give almost the same values. We also checked that vary- 



ing laser beam widths and detunings does not affect the 
final result and that it scales properly with gas pressure. 

We suppose that our method could be incorporated 
into experiments in which diffusion is the limiting fac- 
tor, such as GEM, EIT, quantum memories or collective 
Raman scattering, by relatively straightforward modifi- 
cations such as adding a pump laser or a far field camera. 

We have measured diffusion coefficients of Rubidium 
in Neon, Krypton and Xenon. Reliable values of diffu- 
sion coefficients in these gases facilitate setup design and 
data interpretation in the multimode quantum storage 
experiments. Moreover we recommend the use of Xenon 
as a buffer gas in case of Raman interaction and, to the 
best of our knowledge, we provide the first experimen- 
tal value of diffusion coefficient in this gas. We believe 
that this value can also help develop experiments with 
hyperpolarized Xenon. 
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